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Markovianizing Cost of Tripartite Quantum States 

Eyuri Wakakuwa, Akihito Soeda and Mio Murao 


Abstract —We introduce and analyze a task that we call Marko- 
vianization, in which a tripartite quantum state is transformed to 
a quantum Markov chain hy a randomizing operation on one of 
the three subsystems. We consider cases where the initial state 
is the tensor product of n copies of a tripartite state and 

is transformed to a quantum Markov chain conditioned by 
with a small error, using a random unitary operation on A". 
In an asymptotic limit of infinite copies and vanishingly small 
error, we analyze the Markovianizing cost, that is, the minimum 
cost of randomness per copy required for Markovianization. 
For tripartite pure states, we derive a single-letter formula for 
the Markovianizing costs. Counterintuitively, the Markovianizing 
cost is not a continuous function of states, and can be arbitrarily 
large even if the state is close to a quantum Markov chain. Our 
results have an application in analyzing the cost of resources 
for simulating a bipartite unitary gate by local operations and 
classical communication. 


I. Introduction 

Tripartite quantum states for which the quantum conditional 
mutual information (QCMI) is zero are called short quantum 
Markov chains, or Markov states for short ||T|. They play 
important roles, e.g., in analyzing the cost of quantum state re¬ 
distribution a & investigating effects of the initial system- 
environment correlation on the dynamics of quantum states 
0, and computing the free energy of quantum many-body 
systems Q. 

In analogy to the quantum mutual information (QMI) of 
a bipartite state quantifying a distance to the closest product 
states, it would be natural to expect a similar relation between 
QCMI of a tripartite state and Markov states. However, this 
conjecture has been falsihed |j^(see also |j^, Q). The recent 
results show that the relation between QCMI and Markov 
states is not so straightforward ©-Og, particularly when 
compared to the relation between QMI and product states. 

From an operational point of view, QMI quantihes the 
minimum cost of randomness required for destroying the 
correlation between two quantum systems in an asymptotic 
limit of infinite copies 0- This fact and its variants including 
single-shot cases are called decoupling theorems, and have 
played a significant role in the development of quantum 
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information theory for a decade 0-0. In a simple analogy, 
one may ask the following question; Is QCMI equal to the 
minimum cost of randomness required for transforming a 
tripartite state to a Markov state? 

In this paper, we address this question, and answer in the 
negative. We derive a single-letter formula for the “Marko¬ 
vianizing cost” of pure states, that is, the minimum cost of 
randomness per copy required for Markovianizing tripartite 
pure states in the asymptotic limit of inhnite copies. The ob¬ 
tained formula is not equal to QCMI, or not even a continuous 
function of states. Moreover, the Markovianizing cost of a state 
can be arbitrarily large, regardless of how close the state is to 
a Markov state. In the proof, we improve a random coding 
method using the Haar distributed random unitary ensemble, 
which is widely used in the proof of the decoupling theorems, 
by incorporating the mathematical structure of Markov states. 

There are two ways for dehning the property of tripartite 
quantum states being “approximately Markov”: one by the 
condition that the state is close to a Markov state, on which our 
dehnition of Markovianization in this paper is based; and the 
other by the condition that the state is approximately recover¬ 
able 0, i.e., there exists a quantum operation £ : B ^ BC 
such that « £{p^^). Ref. proved that the latter 

condition has a direct connection with QCMI, namely, small 
QCMI implies recoverability with a small error. 

In fTSl , we introduce another formulation of the Marko¬ 
vianizing cost by employing the concept of recoverability, and 
prove that the cost function is equal to the one obtained in this 
paper for pure states. We then apply the results in analyzing 
the cost of entanglement and classical communication for 
simulating a bipartite unitary gate by local operations and 
classical communication 0. As a consequence, we prove in 
1201 that there is a trade-off relation between the entanglement 
cost and the number of rounds of communication for a two- 
party distributed quantum information processing. 

The structure of this paper is as follows. In Section 
we review mathematical theorems regarding the structure of 
quantum Markov chains, which are extensively used in this 
paper. In Section III we introduce the formal dehnition of 
Markovianization, and describe the main results. Outlines 


of proofs of the main results are presented in Section IV 


In Section |y] we describe properties of the Markovianizing 
cost. In Section [Vl] we calculate the Markovianizing cost 
of particular classes of tripartite pure states to illustrate 
its properties. Conclusions are given in Section VII See 
Appendices for detailed proofs. 


Notations. A Hilbert space associated with a quantum 
system A is denoted by T-L^, and its dimension is denoted 
by dA- For p G SiJ-L^), we denote supp[/9] C by T-L^. A 
system composed of two subsystems A and B is denoted by 
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AB. When M and N are linear operators on T-L^ and , 
respectively, we denote M ® N as ® for clarity. We 
abbreviate ® \4>)^ as |7/))^|(^)®. The identity operator 
on a Hilbert space is denoted by I. We denote [M^ ® 
as and (M^ 0 ® /^)t 

as We abbreviate ® as 

. When f is a quantum operation on A, we 
denote [S ®'vi^){p^^) as {S^ ®id^){p^^) ox S^{p^^). For 
pAB^ pA represents Axb[p'^^\- We denote \'4>){'4>\ simply as '0. 
A system composed of n identical systems of A is denoted 
by A" or A, and the corresponding Hilbert space is denoted 
by or 'H^. The Shannon entropy of a probability 

distribution is denoted as H{{pj}j), and the von Neumann 
entropy of a state p^ is interchangeably denoted by S{p"^) 
and S{A)p. logx represents the base 2 logarithm of x. 


IT Preliminaries 


2) A quantum operation E on S{'Hq) leaves all pk invariant 
if and only if there exists an isometry U : Hq -a- Bq ® 

such that a Stinespring dilation of £ is given by 
£{t) = Txe[UtW], and that U is decomposed by T as 

(T (g) i^)UT^= b')0T" ® ® (2) 

3&J 

Here, are the identity operators on 'Hj'’ := 

suppX^fe P]\k 7 and Uj : 'Hj‘‘ -A 0 B^ are isome¬ 
tries that satisfy TxE[UjUjjUj] = ojj for all j, where 
7^“^ := supp ujj. 

3) r satisfies 

imgT = 0 7f ““ 0 (g -Hf, (3) 

3&J 

where BJ' are one-dimensional subspaces of B°‘'' 
spanned by \j). 

4) B°'“, B°'^ and 77“^ are minimal in the sense that 


In this section, we present a decomposition of a Hilbert 
space called the Koashi-Imoto (KI) decomposition, which is 
introduced in HD and is extensively used in the following 
part of this paper. We then summarize a result in Q, which 
states that the structure of Markov states is characterized by 
the KI decomposition. 


A. Koashi-Imoto Decomposition 

For any set of states on a quantum system, operations on 
that system are classified into two categories: one that do not 
change any state in the set, and the other that changes at least 
one state in the set. It is proved in |2I| that there exists an 
effectively unique way of decomposing a Hilbert space into 
a direct-sum form, in such a way that all quantum operations 
that do not change a given set of states have a simple form 
with respect to the decomposition. We call this decomposition 
of the Hilbert space as the Koashi-Imoto decomposition, or the 
KI decomposition for short. As we verify in Remark in this 
section, the KI decomposition is equivalently be represented in 
the form of a tensor product of three Hilbert spaces. Theorem 
3 in HD^ which proves the existence of the KI decomposition, 
is described in this tensor-product form as follows. 


dim77“" = max dim 77“", dim77“« = max dim 77“" 

ia.T J ia.J J 


Vj € J, 3k s.t. > 0. 

We call F as the KI isometry on system A with respect to ©. 
The KI decomposition and the corresponding KI isometry are 
uniquely determined from ©, up to trivial changes of the basis 
(Lemma 7 in ||2T[). The dimensions of 77“", 77“" and 77“" are 
at most dA- An algorithm for obtaining the KI decomposition 
is proposed in HD- 

It is also proved in HD states {pj\k\k in 

0 are irreducible in the following sense. 

Lemma 2 (Corollary of Lemma 6 in pT| ) The set of states 
{Pj\k^k in 0 satisfies the following properties. 

1) If a linear operator TV on 77“" ;= supp^^fc satisfies 
Po\kNpj\k = Pj\kP 3 \kN for all k, then TV = c/“" for 
a complex number c, where J“" is the identity operator 
on 77“". 

2) If a linear operator TV : 77“" -A 77“/" (j 0= j') satisfies 
P 3 \kNpj\k =Pj'\kPj'\kN for all k, then TV = 0. 


Theorem 1 ( | |2l] , see also Theorem 9 in 0) Consider a 
quantum system A described by a finite dimensional Hilbert 
space B^. Associated to any set of states © := {pk}k 
on A, there exist three Hilbert spaces 77“", 77“", 77“", an 
orthonormal basis {|j)}jGj of 77“" and a linear isometry F 
from 77g := supp(^^ pk) C B^ to 77“" © 77“" (g 77“", such 
that the following three properties hold. (For later convenience, 
we exchange labels L and i? in the original formulation.) 

1) The states in © are decomposed by F as 

FpfcFt = ^ pj\k\j){jr ® ® p% (1) 

3&J 

with some probability distribution {pj\k}jGJ on J := 
{I,-- - ,dim77“"}, states ujj S 5(77“") and pj\k G 
5(77“"). 


Let us now describe an extension of the KI decomposition 
to a bipartite quantum states, which is introduced in |jD- 
Associated to any bipartite state g S{B^ 0 B^'), there 
exists a set of states on A to which system A can be steered 
through ti/AA ^ j g ^ jjjg ^gj states that can be prepared 
by performing a measurement on A' on the state and 

post-selecting one outcome. The KI decomposition of A with 
respect to the set is then associated to . It happens that 
any quantum operation on A which leaves all states in the set 
invariant also leaves invariant, and vice versa. Hence the 
set of operations preserving t^/AA completely characterized 
by the corresponding KI decomposition. More precisely, we 
have the following statements. 

Definition 3 Consider quantum systems A and A described 
by finite dimensional Hilbert spaces 77and B^ , respectively. 
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The KI decomposition of system A with respect to a bipartite 
state g is defined as the KI decomposition 

of A with respect to the following set of states, i.e., 

using C{'H^ ) to denote the set of linear operators on , 

e^A'^A := {p G sin^) 1 3M G cin^') s.t. 

p = TrA'(4) 

The KI isometry on system A with respect to is defined 

as that with respect to &^,ai^a. 

Lemma 4 (See the proof of Theorem 6 in |jT] and Equality 
(14) therein.) Let T be the KI isometry on A with respect to 
, and define := supp[T'^] C LL^. V satisfies the 
following properties. 

1) r gives 

:= ^ Y,P3\j){jr ® ® 

j&J 

(5) 

with some probability distribution {pj}j^j, orthonormal 
basis {|j)}jGj of 71°"', states ujj G and pj G 

2) A quantum operation £ on S{'H^) leaves invariant 

only if there exists an isometry U : —>■ (S> 

such that a Stinespring dilation of £ is given by £^) = 
"Ive\UtU‘^], and that U is decomposed by T as 0, in 
which case we define := suppi/?^". 

We call (|^ as the KI decomposition of on A. The 

following lemma, regarding the equivalence between the KI 
isometries of two bipartite states, immediately follows. 

Lemma 5 The following conditions are equivalent when 

njA _ ijA _ IjA. 

^ ' 

1) A quantum operation on A leaves a state invariant 

if and only if it leaves a state invariant. 

2) The KI isometries on A with respect to and 4'^"^ 
are the same. 

We define the sub-KI isometries as follows. 

Definition 6 Consider a bipartite state , three Hilbert 
spaces and let T be a linear isometry from 

to 7^“" 0 T-L°“- 0 77“®. We call T as a sub-KI isometry on 
system A with respect to if it satisfies Condition 1) in 

Lemma |4l 


B. Markov States 

A tripartite quantum state is called a Markov state 

conditioned by B if it satisfies I{A : C|i?)T = 0. It is proved 
in IJj that the structure of Markov states is characterized by 
the KI decomposition as follows. 

Theorem 7 (See Theorem 6 in |[T| and the proof thereof.) The 
following three conditions are equivalent: 

1) is a Markov state conditioned by B. 


boo l*X*l 


bh » bn 



Fig. 1. A graphical representation of the Markov decomposition of the Markov 
state Each vertex con'esponds to a quantum system, and the white circle 
represents a ‘classical’ system, where the state of the whole system is diagonal 
with respect to The dotted lines represent mixed states. The whole state 
is the probabilistic mixture of the above state with probability qi, namely. 


2) There exist three Hilbert spaces 77^", 77^®, 77*'® and a 

linear isometry T from 77f := supp[T®] to 77*'" 077*'® 0 
77*'® such that is decomposed by T as 

pB^ABCpts = q^\i){i\^“ 0 0 (/)*’®^ (6) 

i 

with some probability distribution {qi}i, orthonormal 
basis {|z)}i of 77*'", states ai G S{H^ 0 77*'®) and 
fii G 5(77*'® 0 77**^). 

3) is decomposed in the form of ® with L being 
the KI isometry on B with respect to 

4) There exist quantum operations TZ from B to BC and 
TZ' from B to AB such that 

T^bc ^ = 7^'(T^'^). (7) 

We call (j^ as a Markov decomposition of a Markov state 
TABC (Figure [^. 

Remark: The KI decomposition is first proved in m by 
an algorithmic construction, and by an algebraic proof in |fT| 
afterward. A similar decomposition is derived in p2] and ||^ | 
in the context of “information preserving structure”. In these 
literatures, the decomposition is given in the form of the direct 
sum of Hilbert spaces as 0^ 77^ 0 77j^. This is equivalent 
to the decomposition in the form of a tensor product of three 
Hilbert spaces described in this section, as verified by choosing 
77“", 77“® and 77“® such that dim77“" = \J\, dim77“® = 
maxj^j'hf and dim 77“® = The corresponding 

KI isometry is defined as 

r:=^|j)“"0(r,.,,0r^0p, (8) 

where Tlj ■ Tif —>■ 77“® and Ti^^j : Hf —>■ 77“® are linear 
isometries and Pj is the projection onto 0 77j^ G 77. 
As stressed in 77“" in Q holds the “classical” part of 
information possessed by pk, 77“® the “quantum” part, and 
77“® the redundant part. 
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III. Definitions and Main Results 
In this section, we introduce the formal definition of Marko- 
vianization, and state the main results on the Markovianizing 
cost of tripartite pure states. The outlines of proofs are given 
in Section Rigorous proofs will be given in Appendix [B] 
andICl 


Definition 8 A tripartite state is Markovianized with 

the randomness cost R on A, conditioned by B, if the 
following statement holds. That is, for any e > 0, there exists 
Tie such that for any n > n^, we hnd a random unitary 
operation V„ : t i-A 2“"^ ™ ^ Markov 

state conditioned by i?” that satisfy 




< e. 

1 

(9) 

The 

Markovianizing cost of 

pABC jg defined 

as 

Ma\b 

[pABC^ _ 1 pABC 

is Markovianized 

with 


the randomness cost R on A, conditioned by B}. 

The following theorem is the main contribution of this work. 
The outline of the proof is given in the next section. 

Theorem 9 Let ^ p^J.g state, and let 

leJ 

be the KI decomposition of on A. Then we have 

j&J 

Based on this theorem, it is possible to compute the Marko¬ 
vianizing cost of pure states once we obtain the KI decom¬ 
position of its bipartite reduced density matrix. However, the 
algorithm for obtaining the KI decomposition, which is pro¬ 
posed in pT[ , involves repeated application of decompositions 
of the Hilbert space into subspaces, and is difficult to execute 
in general. 

Below we propose an algorithm by which we can compute 
the Markovianizing cost for a particular class of pure states, 
without obtaining an explicit form of the KI decomposition. 
The algorithm is based on the following theorem, which 
connects the Markovianizing cost of a pure state and the 
Petz recovery map corresponding to the state. Here, the Petz 
recovery map of a tripartite state from A to AC, an 

idea first introduced in |jT), is defined by 

(Vr e 

A proof of the theorem will be given in Appendix [C] 


Theorem 10 Let be a pure state, such that a CPTP 

map E on S{TL^) defined by 

E (11) 


is self-adjoint. Define another CPTP map Eao by 


1 ^ 

E^-.= lim — Vf", 
TV-s-oo N ^ 


( 12 ) 


n—1 


b) • bo 



Fig. 2. A graphical representation of the KI decomposition of tri¬ 
partite pure states Gg Each vertex corresponds to a quantum system. 
The solid lines express pure states. The whole state is the superposi¬ 
tion of the above states with the probability amplitude y/pj, namely. 


and consider the state 



(13) 

Then we have 



(14) 


Due to this theorem, the Markovianizing cost of pure 
states can be computed by the following algorithm, based 
on a matrix representation of CPTP maps. Here, {\k)}'l^-^ 
is an orthonormal basis of R^, and [■]ki,mn denotes a matrix 
element in the kl-th row and the mn-th column. (See also 
Remark in Appendix C-B ) 

1) Compute -dimensional square matrices Ai, A 2 and A 
given by 


\-^'\\kl,mn 

\A2\kl,mn 


= (A|('k^)5|m)(n|(vl/^)5|0, 


= Tr 




AC\ 


and A = A 2 Aj”^, where the superscript —1 denotes the 
generalized inverse. 

2) Check the hermiticity of A, which is equivalent to the 
self-adjointness of E. If it is Hermitian, continue to Step 
3. If not, this algorithm is not applicable. 

3) Compute a matrix Aqo corresponding to E^o, which 
is given by the projection onto the eigensubspace of 
A corresponding to the eigenvalue 1. Then compute 

Aqo - AqqA]^. 

4) Compute given by 

= X! [^^]ki,mu\k){l\'^ O \m){n\-^'. 

klmn 


5) Compute the Shannon entropy of the eigenvalues of 
, which is equal to 


IV. Outline of Proofs of the Main Theorems 

In this section, we describe the outline, main concepts and 
technical ingredients for the proofs of Theorem and 
Detailed proofs are given in Appendix and 

We first introduce an adaptation of the KI decomposition to 
tripartite pure states as follows. 
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Fig. 3. A graphical representation of the state transformation from ^ in 
0 to j in jl9| by a random unitary operation given by jl8| 


Lemifia 11 Let be a tripartite pure state and suppose 

that the KI decomposition of qjj ^ given by 

3&J 

There exists a linear isometry L' : "Hf —>■ 'h}'" 0 0 

that decomposes together with T as 

j&j 

(15) 

where and are purihcations of and 

respectively, and (j|/)^'’ = Sjji. Moreover, T' is the 
sub-KI isometry on B with respect to 


distribution {pj}j, and fl?"^ is the projection onto the condi¬ 
tionally typical subspace of conditioned by j. Consider a 


unitary operation on supp 


of the form 


y ■= IjXjT 

jeJr,,S 




(18) 


where is the identity operator on supp and Vj” is a 
unitary on the support of 11^"^. We apply on g by 

independently choosing from the Haar distributed random 
unitary ensemble for each j. By this random unitary operation, 
the state © is transformed to the following state 

jeJ„,s 

(19) 


where p'. := Tr,jn“;',|<p,-)(v3^|] and n]- = n“yTr[n“«,]. 
'l '„^5 is a Markov state conditioned by B (Figure!^. 

To Markovianize ^ jj sufficient that we ap¬ 

proximate the transformation from to ( [T^ by V„ with a 
vanishingly small error, where I 4 in are unitaries which are 
decomposed by T®" as ( [l^ . By a random coding method and 
the operator Chernoff bound pT) , it is shown that a sufficient 
number of unitaries in Vn for this approximation is almost 
equal to the inverse of the minimum nonzero eigenvalue of 
(I^, and is given as H{{pj}j^j) + 2J2j^jPjS{p‘'j^) per 
copy. We note that the error e converges exponentially with n 
to zero. 


Proof: The existence of T' follows from Uhlmann’s theorem 
(p4|, see Appendix A-A I. It is straightforward to verify that 
r' is the sub-KI isometry, since we have 


jeJ 

from ( [T5] |. ■ 

J I ABC 

We call (|15[ ) as the Kl decomposition of y) on A and B 
(Figure 1^, and denote it by I'l'ir/). 


A. For Theorem Achievability 

The direct part of Theorem [9] is formulated by the following 
inequality; 

< H{{p,},^j) -b 2^py(y«). (16) 

J&J 

The outline of the proof is as follows. The state |vl/"4BC^®n 
local unitarily equivalent to I'Fkj)®", which is almost equal 
to the state dehned by 


for sufficiently large n. Here, we have introduced notations 
3 = ji - ■ ■ jn, Pj = iFii C) ■ • • 0 Pj„ and ojj = 0 • • • 0 . 

Jn ,5 is the i5-strongly typical set with respect to the probability 


B. For Theorem ^ Optimality 

The converse part of Theorem is formulated by the 
following inequality; 

> H{{P3}3^j) + 2Y,P3S{pT)- (20) 

3^J 

Let US hrst assume tentatively that a Markov decomposition 
of in ^ is given by 

(p/B)®nYABC(p/B)®n ^ ^ p'j\j){j\^" (g) af^ g 

( 21 ) 


with F' being the KI isometry on B with respect to 
In this case, it is not difficult to show that the amount of 
randomness per copy required for transformin g 
to is bounded below by the R.H.S. of (20 1 . Indeed, 

in order to transform I'Fif/)®" to a Markov state in the 
form of (ID, it is necessary that (i) the off-diagonal terms 
with respect to |j) vanish, and (ii) the correlation between 
dp and bpC in the state is destroyed for each 

j. An optimal way for satisfying these two conditions is 
transforming the state GD close to a state of the form 
©■ Since the entropy of the state is approximately 
equal to n{H{{pj}j^j) + the cost of 

randomness required for this transformation is at least about 

However, it might be possible in general that the amount of 
randomness can be further reduced by appropriately choosing 
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£r^ £i £2 £3 £r* 

X L ^ 



Fig. 4. A graphical representation of the channel Due to the completely 
dephasing channel denoted as DP, the system do has some capacity to transmit 
classical information, but has no capacity to transfer entanglement. The system 
dfi has some capacity to transfer entanglement. 


T^bc jjjg corresponding KI decomposition of B. We shall 
see that our choice presented above is indeed optimal. At the 
core of the proof lies the following lemma. 


Lemina 12 Let ^ bipartite quantum state, and let 

: T-L^ —>• 0 0 be the KI isometry on A with 

respect to For any n and e > 0, let be a state that 

satisfies 




AC 


< e, 


( 22 ) 


and let F^^ : —>■ 'H°‘° 0 0 be a sub-KI isometry 

on A with respect to Denoting the decompositions 

pAy[;,4Cp1^ and by and x^kc respectively, 

we have 


- {S{ao)xsKi + 2-S'(aB|ao)x,K/) 

> 5'(ao)^^j + 2S'(aif|ao)>i'„ — C^(e)log(i^. (23) 


Here, C^(e) is a function of e > 0 and 'h, which does not 
depend on n, and satisfies lime_j.o C^r) = 0- See Equality 
(|90|l in Appendix |B-D| for a rigorous definition. 


Proof Outline for Lemma \I2\ Assume here for simplicity 
that 


='H^= (24) 

and suppose the decompositions of qJ ^ac ^g 

given by 


xfSi = E ® C ® (25) 

i 

respectively. Consider a quantum channel on A defined by 




Ft 



which is decomposed as fpt o £3, o £2 ° £1 ° £t ■ The maps 
£y- and are isometry channels corresponding to F^ and 

X i A 

F^, respectively; £i is discarding of system d^; £2 is the 
completely dephasing channel on do with respect to the basis 
|i); £3 is appending of the state conditioned by do (Figure 





c 


Fig. 5. A graphical representation of the state transformation of under 
Ex- The channel E^ has little effect on the state 1 /;®" on average. In 
particular, it almost conserves the correlation that i/;®" initially has. Thus 
the intermediate state has the same amount of correlation due to the 
monotonicity. 


1^. The linearity and the complete positivity of 8 ^ immediately 
follows from ( |26| l, and the trace-preserving property results 
from the fact that F^ satisfies Condition ([^. 

The state x^^ is invariant under the action of 8 ^, and thus 
is almost unchanged due to (22i. By extending the 


data compression theorem for quantum mixed-state ensembles 
1 25 1 , it follows that any state of the form is almost 

unchanged by 8 ^ on average, as long as holds and 

the KI isometry on A with respect to is equal to F^. 

We consider such that its KI decomposition on A is, up 
to an additional decomposition on C', given by 

= Ep^iXiT" 

j&j 


Vki 




® IjXjT'S 


where is a purification of The correlation be¬ 

tween A and C' in the state measured by QMI, is equal 
to 

nI{A:C')^ = 71 /( 000 ^, 0 ^ : c'4)^^^ 

= n(S'(oo)v.Kr + 25'(ofl|ao)^ofJ . (27) 

It can be shown that this amount of correlation is almost 
conserved under 4 . 

Due to the monotonicity of QMI, it follows that the cor¬ 
relation between doaicijj and C' is approximately equal to 
(27 1 at any intermediate step of 8 ^ (see Figure 1^. After the 
action of the completely dephasing channel £ 2 , the system do 
holds no quantum correlation with other systems, and thus 
the correlation between d^dn and C' is bound to be at most 
5'(do) + 2S'(dfl|do) (Figure |^_. Moreover, the state on dodi; 
after £2 is almost equal to x^Pk? *721. A more detailed 

argument reveals that 

£{^ 3 )xski + 2<S'(dfl|do)xajf7 

and consequently proving (j 2 ^. 

C. For Theorem 1701 

Let us first express ( 0 , ([Tgi and 0 in terms of the 
“decomposed” Hilbert space 7^“" 0 0 77“^. A Kraus 


representation of a map £ defined by (111 is given by 


£{■) = with the Kraus operators 

Eu := 


(28) 
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is given by (lOl. For 


Let r be the KI isometry on A with respect to and 

suppose the KI decomposition of 
each E}^i, we have 

Eu ■■= TEkiT^ = Y, IjXjT 

3&J 






(29) 


where 


(30) 


By an extension of Lemma it follows that {ej^ki}ki is 
irreducible in the sense that it satisfies Property 1) and 2). It is 


straightforward to verify from (30i that maps Sj" on 
dehned by £“«(•) := Efe.; (X = 

are trace-preserving. Representations of £ and £ao in the 
decomposed Hilbert space are given by 

£i-)=YEkii-)El 


(31) 


kl 


and Soo '= liinn-).oo(l/-^) respectively, and that of 

is given by £oo{\'^ ki){'^ki\)- 

Due to £ o£^ = and the irreducibility of {ej^ki}ki, we 
have 


3^J 


) UJj 




(32) 


where {p'}jg j is a probability distribution and Cjj {j G J) are 
states on aibobitj^C. Explicit forms of Uj and p' are obtained 
as follows. First, from ( [2^ , (31 1 and the trace-preserving 
property of £j, we have 

Tv[{jr£if)\jr>] = Tr[(id“^®£;«)((jr“f|jr'’)] 

= Tr[or''fbr”] 

for any j and f = TrL^ (r S S{'H^)). This implies that 
the probability amplitude with respect to the basis is 

conserved by £, as well as by £oo- Thus we have p' = pj. 
Observe from (29 1 that £ and £^ do not affect the system a^, 
which implies 


Wi = T"''" ® I jx^r "»■ 

Hence the von Neumann entropy of which is equal to 

that of (p^, is given by 


bsC 


') + 5(7rj' 


5(vl/^«^) = H{{p,),) +Yp3 


In addition, that the self-adjointness of £ implies •P? = 
Since is a purification of , we finally obtain 

that 


b) • bg 



Fig. 6. A graphical representation of a decomposition of tripailite pure states 
for which the Markovianizing cost is equal to QCML The whole state is the 
superposition of the above states with the probability amplitude namely. 


A. General Properties 

Let be an arbitrary tripartite state on finite dimen¬ 

sional quantum systems A, B and C. The Markovianizing 
cost of p^^'^ satishes 

I{A:C\B)p< (p^^^) < I {A : BC)p. (33) 

The second inequality directly follows from the fact that 
decoupling A from BC is sufficient for converting the state to 
a Markov state and that the cost of randomness for decoupling 
bipartite states is asymptotically given by QMI pT) . The 
hrst inequality is proved in Appendix Consequently, the 
Markovianizing cost is equal to zero only for Markov states. 

The Markovianizing cost satisfies a kind of the data pro¬ 
cessing inequality, namely, that 


under any quantum operation £ on C. This is because any 
random unitary operation Vn on A" satisfying (|^ also satisfies 


K (l£^(p)f")- (£^r’^(T 


A’^B'^C’^ 


< e, 


and ^ Markov state conditioned by B". 

As a consequence, an upper bound on the Markovianizing cost 
of a mixed state is obtained as 


Maib(p'^^^) < 

where is a purihcation of p^^. This is because there 

always exists a quantum operation Ep : C ^ C such that 


B. Pure States 


= H{{p,},) + 2 Yp3S{pT) = MA\B{'i’^^^)- 

3 


V. Properties 

In this section, we describe properties of the Markovianizing 
cost of tripartite quantum states. We first consider arbitrary 
(possibly mixed) states, and then focus on the case of pure 
states. 


Let us now consider pure states, based on the result pre¬ 
sented in Section III First, we see that the Markovianizing 
cost Ma\b of two pure states T'l and T 12 are equal if there 
exist A (0 < A < 1) and € S{TiP) such that 


-f (1 - AjT-f O (34) 

This is because = 4/^ and the KI decompositions of A 
with respect to and are equal, the latter of which 
follows from Lemma Indeed, any quantum operation on A 












which keeps invariant also keeps invariant and vice 


versa, as can be seen by observing from (34i that we have 


vk^^*^ = - 


1 


AC 


1 - A 




Second, we obtain a necessary and sufficient condition for the 
Markovianizing cost of pure states to be equal to QCMI as 
follows (Figure |^. 


A. Discontinuity 

We consider tripartite pure states that are expressed as 

l<fA) = 

+\/|^ E w'‘\i=Avf, 

^ fc,Z=l 


Theorem 13 Let be a pure state, and let 

j&j 

be its KI decomposition on A and B (see Lemma [n). Then 
we have = I{A : C\B)qf if and only if there 

exists an isometry F" : -A such that tpj 

is decomposed as 

where {j\j') = 6j^y. 

Proof: We have 


where d = dirnTf^ = dim'H'^, dimAf^ = d^ + 1, 1/d^ < A < 
1, and <1)^ is a maximally entangled state of Schmidt rank d, 
dehned as 

d 

= (35) 

k=l 

For this state, we have 

I{A:C\B)^^ = 2logd-h{X)-{l-X)log{d^-1), 

(36) 

where h denotes the binary entropy dehned by h{X) = 
—A log A — (1 — A)log(l — A). The state is a Markov state 
if and only if A = 1/d^. The distance to the closest Markov 
state is bounded from above by 


I{A : C\B)^ — I{A : C)\ii — I{aoaiaj{ : 
= /(uo ■ 3" I{0‘lQ'R '■ C'|ao)'i'if/ 

= I{aa : + ■ c) 

3 


as well as 

= H{ao) + 2Y,P3SM^^. 


Since is a classical-quantum state, we have /(oq : 

< H{ao) with equality if and only if {supp[(/?^]}j is 
mutually orthogonal. We also have 2S{ait)ip. > I{aji : C)^-, 
which is saturated if and only if 2S{aii)^p. — I{ctR '■ = 

I {(Id : bD)ip- = 0 (see Inequality (|45|l in Appendix A-Bi. 
Flence we have Ma\b{'^^^^) = /(ATc'li?)^ if and only if 


A,ABC AtABC 
“ ^l/d? 


= 2 
= 2 


^1- |(T'A|4'i/d2) 


d^A- 1 


d2 - 1 ■ 


(37) 


The reduced state on AC is given by 
d?X — 1 




AC 


|4'd)($, 


\AC 


1-A 

d2_l 

= A'|$d)($d|^^ + (1-A')^^® 


where tt is the d-dimensional maximally mixed state and A' = 
{d^X — l)/(d^ — 1). Hence the Markovianizing cost does not 
depend on A' when A' > 0, as we proved in Section V-B 


As directly verihed by considering the case of A' = 1, the 
Markovianizing cost is equal to 2 log d for A' > 0. Taking the 
symmetry of T'a between A and C into account, we obtain 


j&J 

which concludes the proof due to Uhlmann’s theorem (p4), 
see Appendix A-A|i. ■ 


An example of states that satisfy the above conditions is given 
in Section IVl-CI 


VI. Examples 

In this section, we consider examples of pure states to 
illustrate discontinuity and asymmetry of the Markovianizing 
cost. We also give an example of states for which the Marko¬ 
vianizing cost is equal to QCMI. 


M^|s(vk,) = Mc|B('I'A) = |^^°®'^ (38) 

Hence the Markovianizing cost is not a continuous function of 
states. In a particular case where A = 2/d^, the Markovianiz¬ 
ing cost grows logarithmically with respect to the dimension 
of the system, whereas QCMI, as well as the distance to the 
closest Markov state, approaches zero as indicated by ( [36] ) and 
(0. 

We note that “approximately recoverable” if A is 

close to 1/d^, i.e., it satishes Equalities Q approximately if 
A « 1/d^. Indeed, since is a Markov state, there exist 

quantum operations TZ : B ^ BC and TZ! : B ^ AB such 
that 


ABC 

1/d 

























9 


Due to the triangle inequality and the monotonicity of the trace 
distance (see Appendix |A-A|l, we have 


< 


< 


A 


ABC 


- 


ABC 


Ijd? 
AB 




ABC 


Ijd? 






ABC 


- 4- 


ABC 

\jtP 



as well as 


Thus Equality p7|) implies 




ABC 


- T' 


ABC 

1/<P 




ABC 


-^('I'a )||i, 




ABC 


, , d 2 A - 1 
< 4 


d2 - 1 


B. Asymmetry 

We consider tripartite pure states that are expressed as 
I'I'a) = 

where d = dimH*^, dimH'^ = diinH^ = d+l, 0<A<1 
and is a maximally entangled state defined by (35i. The 
reduced state on AC is given by 


= A|$d)(4>dr^ + (1 - A) luxur* » tt 


c 


(39) 


Note that does not have any |0)"^|0)‘^ component. 

Hence the CPTP maps on A defined as 0 and ( [T2| ) are 
given by 

£{t) = Tr[PiT]-7ri+Tr[Por]-|0)(0| 

+ ^f’i'rPo + ^-PotPi 
d d 

and 


^oo(r) = Tr[Pir] • tti + Tr[Por] • |0)(0|, 

respectively, where Pq = |0)(0|, Pi = I — Pq and tti = Pi/d. 
It is straightforward to verify that S is self-adjoint. By applying 
£aa to |T'>,) on A, we obtain that 

= A^5^(g)|0)(0|®®^f 

+(1-A)|0)(0|^®|<i>d)($,|^<^. 

Therefore, due to Theorem [T^ the Markovianizing cost is 
given by 


Ma\b{'^\) = h{X) + 2Alogd. 


On the other hand, from ( [39l l, the Markovianizing costs 
Mc\b of does not depend on A when A > 0 as proved in 
Section IV-BI Thus we have 


Mc\B{'d'\) 


2 log d (A > 0) 
0 (A = 0) 


in the same way as (38 i. Hence the Markovianizing cost is not 
symmetric in A and C, as opposed to QCMI, which satisfies 
I[A-.C\B)= I{C ■. A\B). 


C. States for which the Markovianizing cost coincides QCMI 
We consider states that are expressed as 

k=l 


where Afc > 0 and = 1- These states satisfy 

conditions in Theorem [T^ thus the Markovianizing cost is 
given by 


= I{A : = P({A4fc). 


VH. Conclusions and Discussions 

We have introduced the task of Markovianization, and 
derived a single-letter formula for the minimum cost of 
randomness required for Markovianizing tripartite pure states. 
We have also proposed an algorithm to compute the Marko¬ 
vianizing cost of a class of pure states without obtaining an 
explicit form of the Koashi-Imoto decomposition. We then 
have computed the Markovianizing cost for certain pure states, 
and revealed its discontinuity and asymmetry. Our results 
have an application in analyzing optimal costs of resources 
for simulating a bipartite unitary gate by local operations 
and classical communication HE). Some open questions are 
generalization to mixed states, formulation of a classical ana¬ 
log of Markovianization, in addition to finding an alternative 
formulation of Markovianization for which we obtain QCMI 
as the cost function. 

In | |T8| , we have introduced and analyzed an alternative 
formulation of Markivianization and the Markovianizing cost. 
Instead of requiring Condition (|^, we require that the state 
after a random unitary operation is “approximately recover¬ 
able”, i.e., it satisfies Equalities 0 approximately. Eor pure 
states, we have proved that the Markovianizing cost in that 
case is equal to the one obtained in this paper. 
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Appendix A 

Mathematical Preliminaries 

In this appendix, we summarize frequently used facts and 
technical tools used when studying quantum Shannon theory 


and also in the following appendices. Readers who are familiar 
with the material may skip this section. For the references, see 

e-g- 

A. Trace Distance and Uhlmann’s Theorem 

The trace distance between two quantum states p, cr € S{'H) 
is defined by 

--^Wp- cr\\i = \vY sj[p - of 

In the following, we omit the coefficient 1/2 for simplicity. 
For pure states |^), |(/) G TL, the trace distance takes a simple 
form of 

For p,a,T G we have 

llp-r|ii < \\p- (t\\i + \W 

which is called the triangle inequality. The trace distance is 
monotonically nonincreasing under quantum operations, i.e., 
it satishes 

|jp-o-|li > \\£{p) - £{(j)\\i 

for any linear CPTP map S : S{'H) —> S{TL'). As a particular 
case, the trace distance between two states on a composite 
system is nonincreasing under taking the partial trace, that is, 
for p,a G S{H^ (8) TL^) we have 




1 > ||p"^-cr^||l. 


Consider two states p,a G S{TL^) satisfying ||p—cr||i < e, 
and let |'0p)^^ be purifications of the two states, 

respectively. If ds < ds', there exists an embedding of 
into T-L^ , represented by an isometry from TL^ to T-L^ , such 
that 




< 2v^. 


This relation is referred to as Uhlmann’s theorem ( p4| , see 
also Lemma 2.2 in ||29l). In the case of e = 0, the above 
statement implies that all purihcations are equivalent up to a 
local isometry. 


The gentle measurement lemma (Lemma 9.4.1 in |28| ) states 
that for any p G iS('H), X G C{'H) and e > 0 such that 
f) < X < I and Tr[pX] > 1 — e, we have 

y/Xps/X 


P- 


Tr[pX] 


< 2yfe. 


(40) 


As a corollary, when two bipartite states p G S{'H^ 0 TL^) 
and a G S{'H^ ® ) satisfies \\p^ — a^\\i < e, and Her is 

the projection onto supp[tT^] C we have 


..ab 


Tr[p^n^] 

This is because we have 


< 2^6. 


(41) 


> - a" 


= (7 


[A 


-n,p^n,||^ + Tr[p^n/ 

where Il/r denotes the projection onto the orthogonal comple¬ 
ment of supp[(T"^] C TL^, and thus have 


Tr[p^^n^] = Tr[p^n^] = 1 - Tr[p^n;/] > 1 - e. 


^AtjA 


Att-Li 
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B. Quantum Entropies and Mutual Informations 

The Shannon entropy of a probability distribution {px}x^x 
is defined as 

H{{px}x^x) ■= - log Pa;- 

The von Neumann entropy of a quantum state € SifH.^) 
is defined as 

S[p^) = S[A)p :=-Tr[/logp^]. 

If p^ is a probabilistic mixture of pure states as p^ = 
we have S{p^) < Hi{Pj}j) with equality if 
and only if {\(j)j)}j is mutually orthogonal. The von Neumann 
entropy is monotonically nondecreasing under random unitary 
operations, that is, we have S'(A)p < S'(A)v(p) for any random 
unitary operation V on A. For a bipartite pure state , 

we have 

(42) 

For a bipartite state p G S{'H^ ® 'H^), the quantum con¬ 
ditional entropy and the quantum mutual information (QMI) 
are defined as 


S{A\B)p = S{AB)p - S{B)p, 

I{A : B)p = S{A)p + S{B)p - S{AB)p, 

respectively. The von Neumann entropy satisfies the subaddi¬ 
tivity, expressed as 

S[A)p + S{B)p > S{AB)p, (43) 


which guarantees the nonnegativity 
holds if and only if p^^ = p"^ 

Yvjjjcjj js a purification of cr 

we obtain 




AC 


Applying E 




, and by using (4 


S{C),-S{A),<S{AC),. (44) 

Hence QMI is bounded above as 

I{A : C)p < min{2^(A)p, 2S{C)p}. (45) 

For any p G S{'H^ 0 and quantum operation £ on B, 
we have 


S{A\B)p<S{A\B)s^p), I{A-.B)p>I{A:B)s^py (46) 

Inequalities ( [46] l are called the data processing inequality. 

For a tripartite state p G S{'H^ ® 0 B^), the quantum 

conditional mutual information (QCMI) is defined as 

/(A : C\B)p = S{AB)p + S[BC)p - S{B)p - S{ABC)p. 

QCMI is nonnegative because of the strong subadditivity of 
the von Neumann entropy pO) , which is also equivalent to the 
data processing inequality. QMI and QCMI are related by a 
simple relation as 

/(A : BC)p = I{A : B)p + /(A : C\B)p, 

which is called the chain rule. 

For a class of states called the classical-quantum states, the 
quantum conditional entropy and QCMI take simple forms. 


That is, for states p G S{B^ (S> B^) and a G S{B^ (S> B^ 0 
B^), given as 

pXA ^ Y^pi\i){i\^ (g) pf, 

^XAB ^ ^ ^AB^ 

i 

where {|i)}i is an orthonormal basis of B^, we have 
5(A|X)p = 

i 

I{A-.B\X)^ = Y^qJiA-.B),^. 

i 

QMI of a classical-quantum state takes the form of 

I{X:A)p = S{p^)-^p,S{pf), 

where p^ = YhiPiPt- This quantity is equal to the Holevo 
information pi) , and satisfies 

liX : A)p < ^(X)p = 

with equality if and only if {supp[p^]}i is mutually orthogo¬ 
nal. 


C. Continuity of Quantum Entropies 
Define 

f-a;logx (x < 1/e) 

^0(2;) := < 1 / ^ 1 / ^ ’ 

[i (x>l/e) 

f]{x) = X -f f]o{x) and h{x) := r]o{x) -f 770(1 — x), where e is 
the base of the natural logarithm. For two states p and u in a d- 
dimensional quantum system (d < 00 ) such that ||p—cr||i < e, 
we have 

|S'(p) - S'(cr)| < elogd-f 770(e) < 77(e) logd, ( 47 ) 

which is called the Cannes inequality p^ . It follows that for 
two bipartite states p,a G S{B^®B^) such that ||/9—cr||i < e, 
we have 

\S{A\B)p - S{A\B),\ < 77(e) log (dAdl) ( 48 ) 

and 

I J(A : B)p - I{A : B)„\ < 277 (e) log (d^ds). (49) 


D. Typical Sequences and Subspaces ( 1^, see also 

Appendices in ^14^ for further details.) 

Let be a discrete random variable with finite alphabet X 
and probability distribution p^ = Pr{X = x} where x G X. 
A sequence x = (xi, • • • ,x„) G A" is said to be S-weakly 
typical with respect to {px}x^x if it satisfies 

N 

^-n(H(X)+ 5 ) < ^ 

i=l 
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The set of all J-weakly typical sequences is called the S-weakly 
typical set, and is denoted by Tn,s in the following. Denoting 
by Px, we have 

which implies that 

|r„. 5 | < (50) 


A sequence x = {xi,--- ,Xn) G A’" is called S-strongly 
typical with respect to {px\x£X if it satisfies 


^x\x Px 
n ' 



for all a; S A" and = 0 if = 0. Here, is the 
number of occurrences of the symbol x in the sequence x. The 
set of all ^-strongly typical sequences is called the S-strongly 
typical set, and denoted by T* g in the following. From the 
weak law of large numbers, we have that for any e,S>0 and 
sufficiently large n. 


Pr{(Xi,--- ,X„) ern.a> 1-e, (51) 

Pr{(Xi,... ,X„) 1-e. (52) 


Suppose the spectral decomposition of p S 5(7^) is given 
by P = J2xP^\^)(^\- ^-weakly typical subspace 'Hn,s C 
77®" with respect to p is defined as 


Pin,5 ■= span{|a;i) • • • \xn) G 77®"|(a;i,- • • ,a;„) G Tn,s}, 


where Tn,s is the (5-weakly typical set with respect to Px- 
Similarly, the ^-strongly typical subspace 77* g C 77®" with 
respect to p is defined as 


n^g := span{|a;i) •• • |x„) G 77®"|(xi, • • • ,a;„) G 7;*^}. 

Suppose the Schmidt decomposition of G 77"^ 0 77® 

is given by = Tlx \/^\^)^\^)^■ For ^ny (5 > 0 and 

n, let Pin, 5 and 77* g be the (5-weakly and strongly typical 
subspace of (77"^)®" with respect to = TrB[|^/i)(^/i|"^®], 
and let n„ ,5 and H* g be the projection onto those subspaces, 
respectively. From (50i, we have 

rankn„,5 = dim77„,5 < 

From ( |5T| and ( [52l i, we have 

TAKsi\m\^^)n= E 


Tr[n;^,(IV>)(^l^®)®" 

for any e,(5>0 and sufficiently large n. 


xeTn,6 

= E ^ “ r 

.5 


(53) 

(54) 


Appendix B 
Proof of Theorem|9] 

In this Appendix, we show a detailed proof of Theorem 
In the following, we informally denote the composite systems 
aottian by A and bobibjj by B, when there is no fear of 
confusion. 


A. Proof of Achievability (Inequality \16\ ) 

Fix arbitrary n and 5 G (0,1]. Let 5 C J" be the 5- 
strongly typical set with respect to {pj}j^j. For each j G J 
and j = ji--- j„ e Jn,s, define := {l\ji =j,l<l< n}. 
The number of elements in the set is bounded as 

" ■ m) - m) ■ 

For each j G Jn,s, we sort = %°‘«i 0 • • • 0 77“"" as 


(77“")®" = (^ (^ 77“"< 

jtzJ 

For each j and j, let PLj,j,s be the (5-weakly typical subspace 
of PL°'''i with respect to (p“", H^ ^,,5 be the projection 

onto Pij,j,s, and let 11^"^ := J Define 

nf:= E bXjf“ ® ® n“", (55) 

and 

= E (56) 

where we introduced notations pj = x • • • x pj^, pj = 
Pji 0 • • • 0 Pj^ and ujj = ujj^ 0 ■ • ■ 0 . 

Let Vj be any unitary acting on ^j^jPij,j, and define a 
unitary on by 

E b)(jf" ® ® -x (57) 

as 0. We have 

= E (58) 

jeJ„,5 



where \p'j) ■■= Af\pj). 

Let {p{dV),V} be the ensemble of unitaries generated by 
choosing Vj randomly and independently according to the 
Haar measure for each j in •Hzl. Due to Schur’s lemma, as 
an ensemble average we have 


E 




= v 


)p 


.'bsC 


where 7r“" = n“"^/Trn^"^, and 


E 


*7.-5’ 


= 0 


for j f j'. Thus the average state of (581 is given by 


:= E 


K,5(^))('b' 5(^)1 


ABC 


= E Pjbb'XjbT"''"® |wj)(a;j 

jeJ„,5 


10x0 , 


’V 


‘F. 


/6flC 






(59) 
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which is a subnormalized Markov state conditioned by B 
corresponding to (see Figure [^. 

The minimum nonzero eigenvalue of ^ is calculated as 
follows. First, due to the dehnition of J„ 5 , we have 

p. > -^p<P3+s/\J\) ^ 

jeJ 


where 


Suppose Vi, - ■ ■ , Vn are unitaries that are randomly and 
independently chosen from the ensemble {p{dV), V}. Due to 
the operator Chernoff bound (Lemma 3 in pT[), we have 


I ]v ^ i [(1 ~ (1 + eJ'E'n,^] I 


< exp ( — 






Second, since the spectrums of and are the same, 

the minimum nonzero eigenvalue pj of is bounded from 
below as 


Pj 


> 

n 

2-£,.i(S(PU+5) 


jeJ 


> 

n 

2-n(pj+S/\J\)(S{ip°. 


i6J 


> 

2-n 

•(Ej PjS(v“")+i5 log 


where the last line follows from 


E(r + ^) W)+'>) 


- +^(2+log^A) 

3 

3 

Third, we have 

rankn ,-,-5 < 2^^3j\iS{vr)+s) < 2 »(p.+VIJ|)(S««)+5) 

and 


for any G ( 0 , 1 ], which implies that 


Pr 




< 2 e, 


> 1 - 2df exp - 


2”'«Ae2\ 

^ ) 

for an arbitrary R > 0. Therefore, if R satishes 

R > H{{pj}j) + 2j2PjS{Pj'‘) 

3 

+6 {H'{{pj}j) + 21 og (4d^)), 


(60) 


( 61 ) 


and if n is sufficiently large so that the R.H.S. in (j 6 ^ is greater 
than 0 , there exists a set of unitaries such that 


Onfl 


<2e,. 


(62) 


Using unitaries in the set, construct a random unitary operation 
V„ on dndLdn as V„(-) = 2 “"^ Vk{-)V^■ 

Let us evaluate the total error. First, from (j^, @ 

and ( [56] l, we have 

■■= Trpf |vk|7)(vl/|7|] = (vk:, > 1 - (63) 

for any e, > 0 and sufficiently large n. Thus, by the gentle 
measurement lemma (|40ll, we have 


I n,d/\ n,d I 




— 


rankn®;'^ = ]Jrankpj,5 < ]J 
jeJ j&J 

Thus the nonzero eigenvalue Vj of tt?" is, in the same way as 
Pj, bounded from below as 

j&J 

All in all, the minimum nonzero eigenvalue A of ^ is 
bounded as 

A = PjPjVj 

We also have 


rank < \Jn,s\ x rankTr?^ x rank < d^". 


which leads to 


Vn 


-v„(p|7)(vk|7i) 


ddn,S 

Second, from ([62li and (| 6 ^, we have 


^ 2 pp. 


Vn 


I n,d 


){% 


n,S\ 


a 


n,S 




n,S 


a 


n,5 


< 


2e, 


D. 


< 


2e, 


n,S 


1 - e. 


Therefore, by the triangle inequality, we obtain 


v„ (p|7)(3k|7i)- 




Dr, 


^ + 


2U 

1 - e. 


(64) 


From (58 1 and (59 1 , we have Tr[T'„ 5 ] = Dn, 5 , which implies 
that ^n,s/Dn,s is a normalized Markov state conditioned by 
B”. Since the relation ([^ holds for any > 0, i? > 


H{{pj}j) + any 6 G ( 0 , 1 ] that satishes 

and sufficiently large n, we obtain 
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B. Convergence Speed of the Error 

We prove that, in the direct part of Theorem the error 
e vanishes exponentially in the asymptotic limit of n —> c». 
More precisely, we prove the following theorem. 

Theorem 14 There exists a constant c*, > 0 such that for 
any R > sufficiently small 5 > 0 and any 

sufficiently large n, we find a random unitary operation Vn : 
T i-> 2“”^ Yl'k=i on and a Markov state 

conditioned by i?" that satisfy 


C. Proof of Optimality (Inequality ( [2Q| )j 

We assume, without loss of generality, that > dsdc- 
This condition is always satisfied by associating a sufficiently 
large Hilbert space to system A. 

Take an arbitrary R > By definition, for 

any e G (0,1] and sufficiently large n, there exist a random 
unitary operation V„ : t i-G 2“”^ Y^k=i on A and a 

Markov state conditioned by B such that 




< e. 


( 66 ) 


Vn(p®") - 


< 6 exp 




Proof: Let Xi,--- ,Xn be a sequence of i.i.d. random 
variables obeying a probability distribution {px]x- It is proved 
in p5) that there exists a constant c > 0, which depends on 
{Px}x, such that for any d > 0 and n, we have 

Pr{(Xi, • • • , Xn) G Tn,s} > 1 - exp {-cS'^n), 
Pr{(Xi,--- ,Xn) G T*^s} > 1 - exp(-C(5^n). 

As a consequence, there exists a constant > 0 such that 
we have 


Tr[n„,5(|V’)(^/’|^^)®”] = ^ Px > 1 - exp(-c^d^n) 

X&Tn,S 

TrK Px > 1 - exp i-c^S'^n) 


for any d > 0 and n, corresponding to and (54i. Thus, 
for any S > 0, n and E>n,s defined by (|63|, we obtain 

-Dn.5 > 1 - exp (-c^S^n), 


where c,j, > 0 is a constant. Hence we have 


By tracing out A", we have 


< e. 


(67) 


Due to Uhlmann’s theorem (p4), see Appendix A-A I, there 
exists a purification of such that we have 

l^^0n^ABC _ ^ABC < 2y^. ( 68 ) 


Let f' : -A R^" ®R^^ ®R^” be the KI isometry on 

B with respect to ■ From Lemma [TT| there exists a sub- 
KI isometry f : R^ -A R°‘" 0 R°'‘- 0 R°‘” such that the KI 
decomposition of |x) on H and A is given by 

|X„) := (f^®f'^)|x) 

i 

From Theorem]^ and x^^ = a Markov decomposition 

of is obtained by f' as 


= J2q,\i){if“ 0 0 . (70) 


Due to ( 681 and the monotonicity of the trace distance, we 
have 


Vn{\n7){'^%7\)- 


^ O 1 

< 2 exp--— 


Dn,5 


2e. 


1 — exp (—c^(5^n) 






< 6 exp ( — 


,(5^r 




< 2yOe, 


Thus from ([66| and the triangle inequality, we obtain 


for any 6,e, > 0 and n, corresponding to (|^. Substituting 
exp (—c^<5^n/2) into ei, we obtain 


Vn(x^^^) - T 


ABC 


^ ‘2'\/c + e < Sy^. 


(65) 


for any <5 > 0 and n > {hi2)/. 

For an arbitrary i?, choose sufficiently small (5 > 0 such 
that 

R > iT({pj}j) + 2^pjS'((p“«) 

3 

+6 {H'{{pj}j) + 21og (4d^) + . 


Applying fyields 

Vn ^ J ~ ^ Mk 

due to Hence we obtain from 


- T 


ABC 


Vr 


(r^-^x 


ABCPA 


F^ -T 


ABC 

Mk 


<3y/i, 


that 


< 3v^. 


(71) 


Let be the completely depha sing operation on bg with re¬ 
spect to the basis From (70l, we have = 

T)^fV. Thus we obtain from (71 1 that 


Inequality (651 then holds for sufficiently large n, while keep- 


iZ ® . < syi. 


(72) 


ing the R.H.S. in ( [60| strictly greater than 0. This completes 
the proof. ■ 


Here, we defined a random isometry operation 7^ := 
where fpt is an isometry operation corresponding to FV 
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Due to ( |69| l, we have 

2?'''(|X«)(X«I) 

i 

which leads to 

T%c [25'''(IX„)(X„I)' 




Hence we have 




(t; 0 p^'Odx^Xx.. 1)1 = E ® ’ 


where we define 

From ([Toll we have 




y-ABC 
^ Mk 




;\/„'l*'o ^ _A6i 


of being a pure state; the fifth line by Inequality (|74|l; 


where C^'J^e) is a function defined by (901 in Appendix B-D 
Putting together ( |75] l and ( [7^ , we obtain 

— S{ABC)x a S{ao)-irj^j + 2S{a[j\ao)-iii^j 

n 

- ( 27 ?( 3 ye) + C'^(2yfe)) logi^dAdsdc)- 
Noting that is a mixture of 2"^ (not necessarily 


orthogonal) pure states, from ( 661 , we finally obtain 
1 


R > 




(73) 


> -S{ABC)r - vie) log(d^dBdc) 

n 

> S'(ao)\i>^j + 25'(afl|ao)>p^j 

- (3 t ?(3 v ^) + Q2^/e)) logidAdedc) 

= i7(fe},ej) + 2Ep.^(¥’?) 

j^J 

- (3p(3Ve) + C(2ye)) logidAdsdc), 

which implies ( [20| by taking the limit of e —0. 

D. Proof of Lemma \12\ 


Therefore, by tracing out huC in ( |72| l, we obtain 

E *)(d® - E 7*1 *)(d ® 

i i 1 

Thus, by Inequality ( |47| ), we have 

H{{qA^)PY.^^SiafA 

i 

> H{{qi}^) - Vii'fe) log(d^(iB). 

i 

Since the von Neumann entropy is nondecreasing under ran¬ 
dom unitary operations, we have Siffiy,) > SiflA for each 
j from ( |73| l. Hence we obtain 

E 'A E - vi^'/e) logidAdB)- (74) 

i i 

The von Neumann entropy of the state is then 

bounded below as 

SiABC)r = S{AkkbnC)rMu 
= 5(So)t„. + S{AkKC\k)r^, 

= 77({9*}*) + E <l^iSi4"A + Si4>\’>^)) 

i 

= Hi{q^}i) + ^q^{S{af''A + Sif^l^)) 

i 

> H{{qAi) + 2E7i'S'(X’r) “ eivi^y/l) logidAdsdc) 

i 

= SidAxKi + 25'(ajs|ao)x/ir/ “ nviiyfe) logidAdsdc)- (75) 
Here, the third line follows from ([TOli; the fourth line because 


The key idea for the proof of Lemma 12 is similar to the 


one used in ]25) . Let ip^'^ be a state such that the KI isometry 
on A with respect to is the same as that with respect to 
and that it is decomposed as 

:= 






P3\i){jr ® ® ® I7)(jT 


(77) 


where :% 


C 


A- T-p« is an isometry and is 

a purification of (/j®". The state satisfies ip^ = Note that 
we have 

dc = 


3 


rank fj" = 2_^ rank < dA- 

3 


Let € be the set of all linear CPTP maps on SiPL^), and 
define two functions f,g'. £ —M by 

8iip^^')-ip^^' 


fin = 

9i£) = 


Since the Kl decomposition of A with respect to and 

that with respect to ip^'^ ^re the same, f{E) = 0 if and only 
\f g{£) = Q (see Lemma 1^. Define 


C(e) := sup{p(£)|/(£) < e}. 
eee 


(78) 


and the sixth line from ( |69| l. From Lemma ( |68) l implies 

^i^o)xKi + 2S'(aij|ao) Xki 

> n {S{ao)^j^j + 2Sian\aA<i-Ki “ C(2x/e) log d^) ,(76) 


This is a monotonically nondecreasing function of e by defini¬ 
tion, and satisfies lime_>o C<t(c) = 0 as we prove in Appendix 

EE] 

We consider a general situation in which the relation ( |24| ) 
does not necessarily hold. Let H^ be the projection onto 
PL^, and H^ be that onto its orthogonal complement. Using 


a quantum channel on SiPL^) defined by (261, construct 
another quantum channel on S{PL^) by 

£*iT) = £x{^x^^x) + nEiix (Vt g Sin^)). (79) 
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Define quantum channels Si on Ai (1 < / < n) by 

Si{t^') = [f* (g) • • • (g) (g) 

0qf^i+i 0 ••.(g)^'^")] , 


£’*(x"^‘") = £^{x^^) = thus from (22i and 

the triangle inequality, we have 


^■qj0n'^AC _ 


n\AC 


< 2e, 


1=1 


We also have 

n 

= S{C%^. - J2 S{Cl\A^ ■ ■ ■ A^C[ ■ ■ ■ C'_i)£.(v,®0 


/ = 1 


>Y,S{C[)^-Y.S{C[\Ai)s.^^^.) 

n n 

n 


1=1 


1=1 


ill. 


AC' 


Tr[nx(^/>^)®”] 


and 


where Trx\A, denotes the partial trace over 
Ax'■ ■ Ai_iAij^i-■ ■ An- From (25 i, we have 


as depicted in Figure]^ From Condition (221 and tp^ = 
we have — X^\\i < Thus, due to (79i and 

Inequality ( [4T] ), we have 


< 


< 2V- 


which leads to 


I{A : < I{A : + 4nr7(2-ye) log 

(83) 


which implies 

by taking the partial trace. Thus we have 

for any 1 < Z < n. By Inequality ( |4^ and dc < Aa, it 
follows that 

I{A : C')^ - I{Ai : C'i)si(x,) < 4p(C,j(2e)) log 
and consequently, that 

n 

nI{A : C')^ - '^I{Ai : Q')£i(^) < 4n?7(C,,(2e)) logd^- 


by Inequality ( [49| ) and dc < dyt- By the data processing 

(84) 


inequality, we also have 

I{A : < /(oodij : 

Consequently, we obtain from ( |8^ , ( |8^ and ( |84j l that 

nI{A:C% < I{dodn:C%^ 

+471 {r]{2^/e)+ 7]{CA‘^€))) logdA- (85) 

The QMIs in ( [85] l are calculated as follows. First, from ( 
we have 

{S 2 oSio Srjix^) = xTki = E 7*l*X*r'’ ® c 


( 86 ) 


(80) 


Therefore, from the monotonicity of the trace distance under 
£"2 o fi o Sr^, Equality (221 and = 4/'^, we obtain 


^j,^aaR _ apOfl 

+ 71 A, 


< 


{ip' 


^0n\A _ ^A 


< e. 


(87) 


( 81 ) 


Here, we used the fact that £* on A does not change the 
reduced state on C", and that 

Trx\A,.C'\C,' K (X'®”)] = 

because of = ipp. Combining (8O1 and (811, we obtain 


nI{A : C')^ < I{A : + 4np(C(2e)) logd^. (82) 

Define 


^£,oSxo£r)iiP^^') 


VsKI 

Due to ([43, ([481 and 

dim7f“",dim7f“« < dimTf^ = d% 

( |87| ) implies that 

l'S'(ao)x,if/ - >5'(ao)v-,^l < «'7(e)logdA, 
|<S'(dfl|ao)x,„ - S{aR\do)^J < 3np(e) log d^, 

and consequently, that 

S{cio)tp^ + 2iS'(dfl|do)'!/>n 

< S^i^o)xsKi + ‘^S{dii\do)xsKi + 7rir]{e) logd^. 

Since ip^>°-«^ is a classical-quantum state between dp and 
ditC', we obtain 

I{di)d[i : C')^^ 

= i{dc-C%^+i{dn-.c'\dA^^ 

< S{df))'tp„ + 2S'(di;|do)^„ 

< ‘S'(do)x,K/ + 25'(ajf|dn)x,K/ + 7n77(e) log d^ 

= H{{qi}i) + 2'^qiSi(p^’') + 7nr]{e)logdA, (88) 

i 

where the last equality follows from ( [86| ). 

It is straightforward to obtain from (|77[) that 

I{^A . C I{^ciq(1r(Xr . 

= H{{p,},) + 2j2PjSi+';^). (89) 
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Combining ( [85) >, ( [88) > ( [89) >, we obtain 

< H{{q,},) + 2 '^ q,S{(l}i’’) + nC^{e) log dA, 


where 


C(e) = 11?7(2v^) + 477(C^(2e)) 


(90) 


and Csp is a function defined by ( [78| . Thus we finally arrive 
at 


E. Convergence of C,^ 

We prove that C<i/(e) defined by (781 satisfies lime_j.o C<i>(c) = 
0, based on an idea used in | |25l . Due to the Choi-Jamiolkowski 
isomorphism, € can be identified with S{H^ 0 T-L^)- Hence 
£ is compact, which implies that the supremum in ( [78] l can 
actually be the maximum: 

C(e) = max{5(£’)|/(£’) < e}. 


Hence we have that 


Ve > 0, e £ : g{8) = CJe), f{£) < e. 

Define a := limc_>o C«,(e)- Due to the monotonicity, we have 
C«-(£) ^ ct for all e > 0. Consequently, we have that 

Ve > 0, S € : g{£) > a, f{£) < e. (91) 

Define £„ := S C: | g{£) > a}. Due to the continuity of 
g, <£a is a closed subset of <£. Hence 

,3:= mm/(£) 

exists due to the continuity of /. By definition, we have that 

V£ e € : g{£) >a^ /{£) > /3. (92) 

Suppose now that a > 0. We have f{£) > 0 for all £ £ €a 
due to Lemma |5 Thus we have /? > 0, in which case ( |92| ) 
contradicts with (9T]l because e can be arbitrarily small. ■ 


Appendix C 
Proof of Theorem[To1 


In this Appendix, we prove Theorem 10 based on irre- 
ducibility of the KI decomposition. 


A. Irreducibility of the KI decomposition 

Similarly to the irreducibility of the KI decomposition of a 
set of states presented in Lemma the KI decomposition of 
a bipartite state defined by Definition also has a property of 
irreducibility as follows. 

Lemma 15 Suppose the KI decomposition of on A is 
given by 

3&J 


and define Pj^ki ■= {k\^, where {|A:)}fe is an or¬ 
thonormal basis of ■ Then the following two properties 
hold. 

1) If a linear operator N on Hj” := supp[:^“''] satisfies 
PjNpj ki = Pjipj^kiN for all k and I, then N = c/“'' for 
a complex number c, where IJ" is the identity operator 
on -HJL 

2) If a linear operator N : —)■ {j f=- j') satisfies 

PjNpj^ki = Pj>Pj',kiN for all k and I, then N = 0. 


Proof: Define 

PM := 

for M G The set of steerable states corresponding 

to (|^ is given by := m }Hence the KI 

isometry T on A with respect to equal to that with 

respect to by Definition Thus T'm is decomposed 

as 

jeJ 

where V}mg£(w‘^) i^ ^ set of states which is irreducible 
in the sense of Lemma |2] 

To prove Property 1), suppose that N G €{71 j") satisfies 
PjNpj^ki = PjPjM^ foi' ^ii ^ Since is de¬ 
composed as Pj,M = J2k follows that 

PjNpyM = PjPj,MN for all M G €{7-L^). Hence we obtain 
Property 1) due to the irreducibility of 
Property 2) is proved in a similar vein. ■ 


B. Proof of Theorem 10 


Let us first adduce a useful lemma regarding fixed points 
of the adjoint map of a linear CPTP map. 


Lemma 16 (See Lemma 11 in ||T].) Let be a linear CPTP 
map on S{7i), the Kraus representation of which is given by 
£{■) = J2k Dot £* be the adjoint map of £ defined 

by £*{■) = J2k Ell)Ek- Then X e €{71) satisfies £*{X) = 
X if and only if [Ek,X] = [eI,X]=0 for all k. 

The proof of Theorem [^proceeds as follows. Let L be the 
KI isometry on A with respect to and let 

^AC ^ J2Pj\j){jr ® 

3&J 

be the KI decomposition of on A. We have 

3&J 

3GJ Vf'j 

and 

3&J 
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Hence the Kraus operators of E defined by 
as 


It follows from (111 that Due to Lemma 16 

we have 


is decomposed Due to ( |29] l and ( |9^ , this is equivalent to 


3/\'r3'\ 


yk,l; [Eki,^^] = 0 , 

or equivalently, have 

yk,l; [EkiA^]=0, 

from which it follows that 


yj,k,l; = 0. 


Using (301, we obtain 

yj,k,l; (/?“"] = 0, 


where tpj^ki ■= 




\l)^■ Therefore, due to the irre- 


ducibility of the KI decomposition, we have that 
= 7r;« := /;7(dim-H“^), 
and consequently, that 

ej.fc/ = (dimH7)®i^j.fei. 

which implies the irreducibility of {ej^ki}k,i- 

From (281, for any f = FtF'I' (r S 5('H"^)), we have 

(jr»£(f)!/)“'' = EjpmrTijT'). 

where Ejji is a linear map on dehned by 




kl 


Thus we have 


Hence we have 





(J=/) 

(jV/) 


due to the irreducibility of {ej^ki]k,i- 


From (931, we obtain 


= y] Pj I j, J, j, w, 1“"'"“^'’^ 0 

jeJ 


Since £oo(|^if/)(^'tf/|) and are equivalent up to local 

isometries on A and B, we finally obtain 

= Hi{p,}) + Y,Pj 

3 

= Hi{p,}) + 2j2pAAT) 

3 

where we used the fact that = A” is a 

pure state. ■ 


Remark: From 0 and 0- it is straightforward to verify 
that the statement of Theorem does not depend on a partic¬ 
ular choice of a purihcation of That is, for any purihcation 
of T''^, we have = S{ElA^ ), where 

:= )■ A purihcation of is simply 

obtained by 


ursooiAifr = (id“^ ®cyy)(0'r''^ij')“"), 

where Eoojj' is a linear map on dehned by 


N 




:= lim -YE",. 
n.—^nc) Af • ^ 


-'OOjj' ■— 11111 

n—¥oo i\ 

n—1 

Therefore, from 


j&j 


we obtain 


33' 

= Y VAPfAifA ® 

33' 

(93) 

Consider that we have E o E^o = Eao, and thus have 
it follows that 

yk,l; 


and consequently, that 

yk,l; [E^i^I^^,E^{\dkKi)AKi\)]=0. 


dA 

k^l 

and its matrix representation is given by 

= Y Ai]ki,mn\k){Y \m){n\^ . 

klmn 

The matrix elements of the generalized inverse matrix 
are given by 

[Ar']fc/,,nn=(fc|(vE-^)-5|m)(n|(vI/^)-5|Z). 

In addition, the dimension of the eigensubspace of Aoo cor¬ 
responding to the eigenvalue 1 is at least 1, since we have 
Eoo{I) = I due to the self-adjointness of E. These facts justify 
the algorithm described in Section [III| 


Appendix D 

Proof of Inequality ( [33| ) 

The hrst inequality in ( [3^ is proved as follows. For an 
arbitrary n and e > 0, let V„ : r 2“”-^ YAk=i ^kTV^ be 
a random unitary operation on A", and let ^ be a 
Markov state conditioned by such that 


V„(p®")-T 


A^B’^C" 


< e. 


(94) 
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Let be a purification of and be a quantum 

system with dimension 2”^. Defining an isometry W : > 

EA^ by VL = 1^)^ ® ^ Stinespring dilation of 

Vn is given by V„(t) = Tt:e\WtW^. Then a purihcation 
of := Vn{p'^^) is given by 

For this state, we have 

nR > S{E)^,^ = 

= S'(A”B”C'")p; - 

= S'(A”S”C'")p; - nS{ABC)p, (95) 

where the second line follows from ( |44l l. From ( |94| ), we also 
have 


> - np{e) log idAdedc) 

= 5'(A”B")t + S{B^C^)r - S{B^)r 

-np{e) log (dAdsdc) 

> + S{B^C^)p>^ - S{B'^)p,^ 

-4n?7(e) log (dAdsdc) 

> S{A^B^)p^„ + 

-4n77(e) log {dAdsdc) 

= n(S(AB)p + S{BC)p-S{B)p) 

-4np{e)\og{dAdBdc)- (96) 


Here, the second line follows by Inequality ( |47] l; the third line 
because of T being a Markov state conditioned by S"; the 
fourth line by Inequality ( |47| ); and the hfth line by the von 
Neumann entropy being nondecreasing under random unitary 
operations, in addition to From (95i and 

we obtain 


R> I{A: C\B)p - 4p(e) log (dAdsdc), 


which concludes the proof. 



